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Abstract 

We calculate the magnetic moment caused by the persistent currents in polyg- 
onal carbon nanotube tori by the tight binding model. The polygonal CNT 
tori are formed by introducing heptagonal and pentagonal defects along the 
inner hole and outer fringe, respectively. We found a new type of large param- 
agnetic persistent current caused by the semi-metallic band structures of the 
periodic CNT junction. It is contrasted with the persistent current originated 
from the metallic band structure. 

In spite of diamagnetism of the graphite, this paramagnetism is caused by 
not only the magnetic flux in the inner hole (AB flux) but also the magnetic 
flux directly penetrating the graphite plane (direct flux). This magnetic mo- 
ment is close to that calculated with only Aharonov Bohm effect where the 
direct flux is included in the AB flux. 

I. INTRODUCTION 

In spite of semi-metallic nature of the graphite, the carbon nanotube (CNT) becomes 
metallic and semi-conducting depending on the radius and the helicity of the honeycomb 
lattice. [1] This promising feature in nanotechnology is caused by the quantum size effect; 



the system smaller than the coherent length shows electronic states different from those of 
the macroscopic system. Another example of the quantum size effect can be seen in the 
persistent current of the mesoscopic ring. An ideal clean mesoscopic one-dimensional metal 
ring is diamagnetic when the number of the electron is 4n + 2 with an integer n, while it is 
paramagnetic otherwise. [2] It means that the ring can have opposite sign of the magnetic 
susceptibility compared to the bulk system. We should notice that this positive magnetic 
susceptibility is inversely proportional to the temperature, i.e., it becomes divergent at the 
absolute zero temperature. [2-4] Though it seems similar to the Curie Wise rule, the origin 
is not the spin but the persistent current. As there is no general terminology for this 
paramagnetism, we call it giant orbital paramagnetism (GOP) in the present paper. 

We are interested in the CNT ring discovered recently for the following reasons. [5] As 
for the magnetic susceptibility with respect to the magnetic field perpendicular to the hon- 
eycomb plane, the CNTs and the graphite show similar negative values, [6] whereas the 
persistent current of the CNT ring can show the GOP. [4] Furthermore in order to make the 
radius of the ring R nanometer-sized, the CNT ring will be more suitable than the conven- 
tional mesoscopic ring of metals or GaAs/AlGaAS. [7] Motivations to downsize the persistent 
current ring are as follows; (l)The circulating persistent current / equals evF/{27rR) with 
Fermi velocity vp, and is related to the magnetic moment per ring M as M = IttR^. In 
the rings closed-packed, therefore, the magnetic moment per unit area is proportional to / 
and can be increased by shrinking R. (2) The persistent current occurs on condition that 
the coherent length is larger than 27ri?. Since decreases as the temperature rises, the 
persistent current can be observed at higher temperature as R decreases. 

In other theoretical works on the persistent current of the CNT ring, both the 'discli- 
nation' and the 'direct flux' are not considered, [4,8-10] whereas we consider both of them. 
Shrinking the radius of the CNT ring R increases importance of these two factors. Firstly 
it is explained for the dischnation as follows. Figure 1 illustrates a 'polygonal CNT torus' 
where pentagonal and hcptagonal defects form the corner of the polygonal shape along the 
outer fringe and the inner fringe, respectively. These defects are generally called disclinations 
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and have significant effect on the CNT junction, [11] the CNT cap [12] and the hehcal CNT. 
[13,14] Here we define a 'polygonal' CNT torus as the CNT torus with the disclinations, 
though it can also have a rounded shape when the disclinations are close to each other. [15] 
Only few attempts have so far been made at the persistent current of the polygonal CNT 
torus [16] in contrast to that of the 'circular CNT tori' discussed by many authors. [4,8-10] 
Here we define a 'circular' CNT torus as the CNT torus formed by elastic deformation of 
the straight CNT with no disclination. In the circular CNT tori, the honeycomb lattice 
shrinks along the tube axis with a factor Ri/ {Ri + D) at the inner fringe compared to that 
at the outer fringe. Here Ri and D denote the radius of the inner hole and the diameter of 
the CNT, respectively. Though this axial strain is relaxed to a certain degree by the elastic 
deformation called 'buckhng', [17] it will not be enough to keep the stability if i?i is close to 
D. In Fig.l, on the other hand, this strain is relaxed by decreasing the number of hexagons 
along the inner hole side. Though it is widely known that the impurity reduces the GOP, 
[3,8] the effect of the disclination on the GOP is open to question. [3,4,8-10] Secondly, the 
importance of the direct flux is explained as follows. We consider both the magnetic flux in 
the graphite plane and that in the inner hole while other researchers have investigated only 
the latter flux. We call the former flux 'direct flux' and the latter flux 'AB flux' because 
the former touches the electron directly and the latter induces the Aharonov-Bohm effect. 
Ratio of the AB flux to the direct flux is R\ to -\- DY — R\. This ratio indicates that 
the direct flux cannot be neglected when Ri is comparable to D. Because the graphite is 
diamagnetic under the direct flux, we have to find out whether the direct flux reduces the 
GOP. In order to clarify these open questions, we discuss the infiuence of the disclination 
and the direct fiux on the magnetic moment of the polygonal CNT torus. 

II. TIGHT BINDING MODEL 

When the polygonal CNT torus is composed of only semi-conducting CNTs, it has finite 
HOMO-LUMO gap. As can be known later, it means that the GOP does not appear. 
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Thus our discussion can be limited to the metaUic CNTs. As first examples, we discuss 
the polygonal CNT torus with symmetry Dg/i formed out of the armchair CNTs which are 
metallic irrespective of their radii. We use the tight binding (TB) model with only tt orbital. 
The position of the atom is given by the assumed polygonal shape of the torus. 

Figure 2 illustrates a part of the projection map of a polygonal CNT torus. The torus is 
composed of six unit cells, one of which is represented by the rectangle ABB'A'. The prime ' 
means that points X' and X become identical when actual three dimensional shape is formed. 

— * 

Vector L denotes the circumference of the original armchair CNT. The length and direction 
of the CNT axis per unit cell along the outer fringe (inner fringe) is represented by S {S2). 
The i membered ring disclination is formed at Pj and Qj {i — 5,7) by removing the shaded 
area and sticking the lines P5P7, P7D', DQ7 and Q7Q5 on the lines P5P'7, P'7C', CQ'7 and 
QVQs, respectively. The disclinations approximately keep the sp^ bonds of carbon atoms, 
i.e., the number of nearest neighbors is three for each atom. In spite of the unchanged local 
structure, the phase of the wave function is shifted when the electron circulate around the 
disclination owing to excess or deficiency of the tt/S angle on the projection map. This 
topological effect cannot be represented by the effective potential energy or modification of 
the bond strength. 

Because of the Dg/i symmetry, R^^ and i?77 are parallel to L, while ^2 and S are perpen- 

— * — * _ 

dicular to L, where the definition of Rij is clearly illustrated in Fig. 2. Thus the four integer 
parameters, — \L\/{V3a), ^75 = |^75|/a, ns — \S\/a and ^77 = |^77|/(-\/3a) specify the 
CNT torus with a ~0.25 nm being lattice constant of the graphite. For example, these pa- 
rameters are ul — 6,ns — 5, ^75 = 2 and 7177 = 2 in Fig. 2. To keep the bond lengths almost 
constant, the cross section cut along P5P7Q7Q5P5 has to be rectangular, so the four integers 
have to satisfy the condition 77,77 = ('^l ~ ^75) l'^- We choose the uniform magnetic field B 
parallel to the six-fold rotational axis of the torus. In Fig.2, denotes the component of 

— * 

B perpendicular to the projection map. 

The unit cells are numbered along the CNT axis with ni as shown in Fig.l. Figure 2 
shows labels of atoms in each unit cell, (712, n^), where 712 specifies the zigzag rows along the 
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tube axis. At the inner hole side, n2 — 1 (^2 = a-nd n2 increases (decreases) one by one 
as the zigzag row chmbs up (down) the inner hole wall and approaches outer fringe. In 712 
zigzag row, atoms are numbered along the CNT axis as na = 1, 2, • • • , q{n2). For example, 
g(l) = g(2) = 6, q{3) = 7 in Fig.2. Between the neighboring unit cells, atoms (ni, ^2, 1) and 
(ni, 712, 9(7^2)) connect with atoms (ni — 1,77-2,9(^^2)) and {ni + l,7i2, 1), respectively. With 
this label, matrix elements of Hamiltonian between atoms n and m are represented by 

H{B)n^^^ -tn,m(i^p{iP{n,m)B) . (1) 

Each atom n has a constant hopping integral —tn,m — — ^(~ — 3eV) with the three nearest 
neighbors m and tn,m — with all the other atoms m. The effect of the magnetic fields is 
included by Peierls phases (3{n,m)B, [18] where 

(2) 



P{n, 


m) 


= 7D(n, m)Ci + 7As(n, m)C2 




Ci 


= \/37raV(20o) , 




C2 


= V37r|>S2|V(20o) . 




m) 


'^n2,m2 ('^ni,mi+l'^n3, 1*^7713, 9(712) '^ni,mi — 1*^7713, 1*^713 ,9(712)) 


loin, 


m) 


'^7l2,7n2{/('^2) ((^9(712)711+713, 5(712)7711 +7713 + 1 



— (^9(7l2)7ll+7l3,?(7l2)77ll+77l3-l) + 5'(^2)7Ab(^, "T') } , 

with the magnetic flux quantum 0o = h/e and Kronecker delta 5. Here / and g are linear 

or constant as a function of n2] 



n75 

1^2! — — 0.5 




• \n2\ > 1+^75 +7177 

• 7l77 + 1 < \n2\ < 7l75 + 7177 

• |7l2| < ^77 



(3) 



9(^2) 







(4) 



• • • P2I >^+n75 + n-j-j 

^ fM + 2 EHil7i /(i) • • • \n2\ < ^75 + nr,. 
In eq.(2), 7/j(i; j)Ci and 7as(^; i)C'2 come from the direct flux and the AB flux, respectively. 
Equation (2) satisfies the condition that 0o-BI]j=i /3(rz*^*+^'';7i(*))/(27r) equals the magnetic 
flux surrounded by the closed loop n^^^ nP'^ — > • • • n^^) — n^^^ . 



III. RESULT 



Firstly we explain the origin of the GOP of the circular CNT torus with the dispersion 
relation Ei{k, B) of the straight CNT under the uniform magnetic field B perpendicular to 
the tube axis. Ajiki and Ando showed that the direct flux flattens the dispersion lines near 
the Fermi level. [18] This is schematically shown by the dashed curves Ei{k, B) compared to 
the linear dispersion Ei{k, 0) in Fig. 3. When the lattice distortion is neglected, the energy 
level of the circular CNT torus e/j can be expressed by 

ei,^ = Ei{k^{B),B) (5) 

ki{B)^{j + {BSKB/(t>o)}/R (6) 

with the area of the inner hole S'ab and an integer j. Equation (6) represents the AB effect 
while difference between Ei{kj{B),B) and Ei{kj{B),0) comes from the direct flux. Then 
the total energy U {B) is obtained as 

U{B) = 2j2hr>{ei,)ei, (7) 
i,j 

where /fd is Fermi-Dirac distribution function and factor 2 represents spin degeneracy. 
Zeeman effect is neglected here but will be discussed latter. We consider the case of the 
absolute zero temperature, so /fd is replaced by the summation J2i,j=occ limited to the 
occupied states. 

The magnetic moment per torus M is calculated as M = —dU/dB. In order to see the 
effect of the direct flux, we firstly discuss the magnetic moment induced by only AB flux 

Mab = -5AB/(i?0o) E {dEi{k,0)/dk)\k=k,(o) ■ (8) 

l,j=occ 

Since Ei{k,0) is an even function of k, the energy level at k — %(0) and that at k — 
k-j{0) — —kj{0) are occupied at the same time. Their contribution to Mab usually cancels 
each other, because dEi{k,0)/dk is an odd function of k. When the highest occupied level 
(HOL) comes at K and K' corner points, however, this cancellation does not occur so that 
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the GOP is induced. This is illustrated in Fig. 3 where the energy levels and its change 
caused by the AB effect are indicated by circles and arrows, respectively. Without the 
direct flux, the levels move along the linear dispersion line and only the level with negative 
dE/dk is occupied (closed circle) while that with positive dE/dk becomes vacant (open 
circle). Only the former contributes to the magnetic moment and brings about the giant 
orbital paramagnetism (GOP). The occupation of only the state with a negative dE/dk 
corresponds to the generation of the persistent current. The essential point is that the 
corner point is partially occupied (gray circle). If the corner point were fully occupied, the 
open circle level would be also occupied and it would cancel the energy decrease of the closed 
circle level. With the direct flux, however, the closed circle level is located on the dashed 
curve and shifted upward compared to the case without the direct flux, i.e., the direct flux 
reduces the GOP. 

The magnetic moment for the polygonal torus can be calculated almost the same way 
as the circular torus. In this case, we use the phase a instead of the crystal wave number k. 
Rotation by 7r/3 with respect to the six-fold symmetry axis is equivalent to multiplying the 
wave function by the phase factor exp(iQ;) as seen in Fig.l. Though knR/S corresponds to 
a, it should be noted that R cannot be defined uniquely for the polygonal torus. The energy 
levels Ei{a) are the eigen values of the matrix Hi{B)+exp{ia)H2{B)+exp{—ia)H2{B) where 
Hi and H2 refer to the bonds in the unit cell and those connecting the neighboring unit cells, 
respectively. They are obtained from the Hamiltonian matrix H as ifi(n2, 713; 777.2, ms) = 
i/(ni, 772, 773; 771,7772, ma), and iJ2(«2, '"•3; "^2, ^3) = H{ni,n2,n3;ni + 1,7772,7773) with the 
label of atoms (rii, 712, n^) defined in section II. Because the wave function is invariant under 
27r rotation, a takes discrete values aj — (7r/3)(j + (SjS'ab/^o)) U — integer) including the 
AB flux, BSab = -B3\/3| 5*2^/2. It results in discrete levels Ei{aj,B) of the CNT torus. 

On the other hand, Ei{a, B) with continuous a represents the continuous energy spectra 
of the periodic CNT junction made by connecting the unit cell of the CNT torus. In the 
present case, the band structure Ei{k,0) is classified by three parameters, 715/3,7175/3 and 
{ns — 7i75)/3 as hsted in Table I according to Ref. [13]. As shown in Fig. 2, the CNT torus 
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contains two kinds of CNTs. One is the original armchair CNT with the chiral vector L, 
and the other is characterized by different chiral vector L'. When rif^/S is an integer, latter 
becomes also metallic and we proved analytically that types 1 and 2 are semi-conducting and 
metaUic, respectively. [14] Though we obtained only numerical results for the other types, 
we can say that they tend to be semi-conducting as the axial length of the semi-conducting 
CNT segment \L' x R^il/lL'] becomes longer. [19] 
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TABLES 



TABLE I. 



integer 


7175/3 

non-integer 


integer 


5/3 

non-integer 


integer 


non-integer 








{ns - 

integer 


- 2n75)/3 

non-integer 


type 1 


type 2 


type 3 


type 4-1 


type 4-2 


semi-conductor 


metal 


semi-metal 
semi-conductor 
metal 


semi-metal 
semi-conductor 
metal 


semi-conductor 
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The dispersion curves are shown for type 2 periodic CNT junction (metaUic) in Fig.4 and 
for type 3 periodic CNT junction (semi-metaUic) in Fig.5. The index I of Ei{a, B) is defined 
as Ei{a, B) < E2{a, B) < ■ ■ ■ ,< E2n{(^, B), where 2A'" is the number of atoms in the unit 
cell, i.e., N = 2(nsnL — n^inTj) — n\'j. Thus the HOMO band is denoted by Em{ol, B) while 
-E'Ar+i(Q;, B) corresponds to the LUMO band. The circles, lines and arrows in Fig.4 and Fig.5 
have same meanings with those in Fig. 3. Firstly we consider the magnetic moment under 
only AB flux Mab- It is analogous to eq.(8) as 



with C2 defined in eq.(2). The HOMO and LUMO bands of type 2 is similar to those of 
metaUic CNTs, but the band crossing point is slightly shifted from a — ±27r/3 owing to 
the phase shift at the disclinations. [14] Therefore the HOL is fully occupied and the GOP 
does not occur. Under finite B, however, a2 — 27r/3 + C2B comes at the cross point so 
that Mab shows a discrete change from negative to positive as shown by Fig. 6. Contrary 
to it, HOL can be partially occupied in the CNT torus of type 3 and 4-1, on the condition 
that the corresponding periodic CNT junction is semi- metallic. Figure 5 is an example of 
it where Ej^ at a = is vacant and E'at+i at q; = ±7r/3 are partially occupied instead. 
The AB flux lifts the degeneracy of Hogs and only the lowered level at a = 7r/3 is occupied 
causing the GOP. Then does the direct flux reduce the GOP ? Surprisingly, it enhances the 
GOP. The dashed hues in Fig.5 represent the dispersion relation Ei{a, B) under the finite 
magnetic field B = 0.01(l)Q/a^. In contrast to Fig.3, the dispersion lines is shifted along the 
a axis rather than along the E axis. By this horizontal shift, the effect of the AB shift is 
not canceled but enhanced. 

Here we should notice the competition between the GOP and spin paramagnetism. When 
the Zeeman spht q/jIbB {/jLb — Bohr magneton, gr ~ 2) is larger than the GOP spht \MB\ — 
\BdE/dB\ , both the HOLs with positive dE/dk and negative dE/dk become occupied by 
the same spin so that spin magnetic moment, that is 2g^B pcr torus, appears instead of 
the GOP. Since jiB — ta'^/4'o — 3 eV, a ~ 0.25 nm), the GOP is relevant when M is 




(9) 



l,j=occ 
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larger than 4i/(0o/a^). To search the polygonal CNT torus showing the relevant GOP, the 
magnetic moment M(1.5AB) = -{U{2AB)-U{AB))/AB with AB = 0.5 x IQ-^^o/a^, i.e., 
B = 1.5AB ~ 0.24 Tesla, is calculated for two hundred seven kinds of the tori in the range 
of the parameters 3 < ^5 < 12, 1 < ^75 < 6, 3 < < 12 and 1 < ^77 = (ul — n'j^)/2 < 6. 
Among the two hundred seven tori, the number of the tori with M larger than 4ta^/(pQ is 
twenty one. All of them belong to type 3 or type 4-1 and show the GOP. Especially three 
CNT tori show M larger than 2Oto^/0o- In contrast to it, M is much less than 4to^/0o 
for the tori without the GOP. Since MB is approximately the split between the degenerate 
HOL induced by B, temperature has to be lower than MB. When M = 2Qt/ 0?) and 
B — 1.5 AS, MB ~ 6 K, which is achievable in the experiment. 

To show the distribution of M and the enhancement of M by the direct flux, Mab/M 
is shown in Fig. 7 as a function of M for the tori showing the relevant GOP. For all 
the tori, Mab < M, i.e., the direct flux enhances the GOP. Figure 7 also shows that 
{\S\^ I\S2\'^)Mab/M becomes closer to unity than Mab/M where the factor \S\'^/\S2\'^ rep- 
resents the ratio of the AB flux to the full flux. This suggests that the effect of the 
direct flux can be almost included by the AB effect for which the cross section is not 
3\/3|52|V2 but 3a/3|,S|V2. It is equivalent to replacing Ci and C2 in eq.(2) with Ci = 
and C2 = \/37r|5|7(20o)- In other word, the HOMO and LUMO dispersion curve can be 
approximated by 

Ei{a, B) c^Ei{a+ (7r/3)(S5^/0o), 0) (10) 

with the cross section of the graphite plane Sd — 3-\/3(|5'p — |^2p)/2. In fact, the shift of 
a in eq.(lO) can be seen in Fig. 4 and Fig. 5 where the dashed lines represent the dispersion 
relation Ei{a,B = O.Olcpo/a^). It is also the effect of eq.(lO) that the magnetic moment M 
shows discrete change at lower B than Mab in Fig. 6. 
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IV. SUMMARY AND DISCUSSION 



In some of the polygonal CNT tori, the persistent current leads to large positive magnetic 
susceptibility inversely proportional to the temperature, which we call the giant orbital 
paramagnetism (GOP). The magnetic moment per torus can be larger than Bohr magneton 
in the achievable condition, B ~ 0.24 Tesla and T < 6 K. The necessary condition for the 
GOP is that the corresponding periodic CNT junction is semi- metallic, i.e., type 3 or type 
4-1. The magnetic flux is divided into two parts which are the AB flux in the inner hole and 
the direct flux intersecting the graphite plane. In spite of the diamagnetism of the graphite, 
the direct flux enhances the GOP and its effect can be effectively included in the AB effect, 
where the cross section of the AB flux is not the inner hole area, but the area surrounded 
by the outer fringe, as shown in eq. (10) and Fig.7. 

The GOP of the circular CNT torus is caused by the metaUic band structure of the 
straight CNT, whereas that of the polygonal CNT torus originates from the semi-metallic 
band structures of the periodic CNT junction. This difference comes from change of the 
bond network, i.e., removing the shaded areas and sticking their edges in the projection map 
illustrated by Fig. 2. Compared to it, the elastic deformation does not affect whether n and 
m are bonded or not, so it is irrelevant to the GOP of the polygonal CNT torus. Owing to 
this new kind of the GOP, the polygonal CNT torus is a promising nano-structure to realize 
the large persistent current. 
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FIGURES 
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FIG. 1. (a)Three dimensional shape of the polygonal nanotube torus. The applied uniform 
magnetic field is perpendicular to this page. Labels of unit cells ni are shown, (b) The unit cell of 
the CNT torus. The torus is composed of the six unit cells. 
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FIG. 2. Projection map of the CNT torus. The right rectangle is the unit cell. Magnetic field 
perpendicular to this page is denoted by B±. Labels of atoms (71-2,^3) in the unit cell are shown. 
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FIG. 3. The solid lines and dashed curves represent the dispersion relation without and with the 
direct flux, respectively. When the magnetic field B equals zero, the highest occupied level (HOL) 
of the torus is at the crossing of the solid lines as indicated by the gray circle. It is degenerate and 
partially occupied. Under finite B, the degeneracy is lifted and only the lower level (closed circle) 
is occupied while the higher level becomes vacant (open circle). The arrow indicate the change of 
the levels induced by the AB effect. 
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FIG. 4. The band structure of the periodic CNT junction of type 2 with the parameters 
(n^, n^, n75, nss) = (5,7,3,2). The circles, hues and arrows have the same meanings with Fig. 3. 
For the dashed hnes,-B = 0.01(/)o/a^. Inset is the magnified band structure near the highest occupied 
level (HOL), E7v(-27r/3, 0). 
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FIG. 5. The band structure of the periodic CNT junction of type 3 with the parameters 
(ns,ni:, 7775,7155) = (3,7,1,3). The circles, Unes and arrows have the same meanings with Fig.3. 
For the dashed lines, i? = O.Ol^o/o^- 



19 




20 



FIG. 6. The magnetic moments of the same periodic CNT junction as that in Fig.4, i.e., 
{ns, riL, ?^75, ?^55) = (5, 7, 3, 2), as a function of the magnetic field B. Here M and Mab denote the 
magnetic moment induced by full flux and that induced by only AB flux, respectively. The former 
is calculated by the total energy U{B) as M((j - 0.5)AS) = -{U{jAB) - U{{j - l)AB))/AB 
with AB = 0.5 X 10~^(^o/a^) and j = 1, 2, • • • , 100. For Mab, the total energy is calculated by 
the Hamiltonian where 7^ in eq.(2) is replaced by zero. 
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FIG. 7. The ratio Mab/M and {\S\^ I\S2\'^)Mab /M as a function of M for the twenty one 
CNT tori with M larger than 4io? /(pQ. Here M and Mae denote the magnetic moment induced 
by full flux and that induced by only AB flux, respectively. The applied magnetic field B equals 
0.75 X 10^5 0o/a^- The factor |5p/|52p means the ratio of the full flux to the AB flux. 
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